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When hot quark gluon plasma expands and cools down after an heavy ion collision, 
charge conservation leads to non-trivial correlations between the charge densities 
at different rapidities. If these correlations can be measured, they will provide 
information about dynamical properties of quark gluon plasma. 



a,: 

1. Introduction 

<D ! 

e-j , The view to the quark gluon plasma phase (QGP) in heavy ion collisions is 

obscured by events taking place later on. One way around this is to focus on 
conserved charges, such as electric charge, baryon number or strangeness. 1,2 
rS ' Later evolution can only change these charges locally, through pair creation 

■ or annihilation, but charge fluctuations with long enough wavelengths will 

remain unchanged. 

One possibility is to consider the fluctuation (AQ 2 ) of the charge Q in a 
given rapidity window A-q. 1 ' 2 Because of charge conservation, this quantity 
can only change when charges move in and out of the window, but if the 
window is wide enough, this effect should be small. Therefore, the measured 
charge fluctuation should reflect the initial value and be different for QGP 
and hadronic phases. 

Ignoring long-range Coulomb forces, one can estimate that in the QGP 
phase, where the electric charges of the quarks are either 2/3 or —1/3, the 
fluctuation is (AQ 2 ) « 0.2A c h- In the hadronic phase, elementary charges 
are ±1, and as a consequence the fluctuation is much larger (AQ 2 ) ss 
0.7A c h- 2 In principle, this difference could be used to find out whether QGP 
formed in the early stages of the collision. However, current experimental 
results are consistent with the hadronic value. 3 - 4 ' 5 They are also very close 
to (AQ 2 ) = A c h, which would correspond to a purely random distribution 
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of ±1 charges. 

In this talk, I will further explore the information obtainable from charge 
fluctuations. In particular, I will consider correlations between charge den- 
sities at different rapidity values, and show that they carry quantitative 
information about different stages of the collision. It is hopefully easier 
to subtract the contributions due to phenomena taking place in the later 
stages of the collision, such as decays of hadronic resonances, from these 
correlations than from the charge fluctuation signal (AQ 2 ). 

2. Diffusion in an expanding system 

In the first stages of a heavy ion collision, the expansion in the direction of 
the beam is much faster than that in the orthogonal directions. The effects 
I will be discussing are all due to this expansion, and therefore I will not 
consider the orthogonal directions. I will also assume that the two nuclei 
are moving at the speed of light so that the collision event is boost invariant. 
In that case, it is convenient to use the Bjorken coordinates t and r\ defined 
by t = Tcosh?;, z — rsinh?]. In these coordinates, the Minkowski metric 
becomes ds 2 = dr 2 — r~ 2 dz 2 , which is the metric of a 1+1 dimensional 
FRW universe with scale factor a(r) = 1/r. This means that very similar 
considerations apply to charge considerations in the early universe, as well. 

I will assume that the evolution of the charge density is purely diffusive. 
Charge annihilations in that case have been studied for a long time, 6 but 
often without coupling the system to a heat bath. In the presence of a ther- 
mal bath, the diffusion equation in Bjorken coordinates for the comoving 
charge density p = dQ/drj is 

d T p=£^d 2 ~p + d n Z v , (1) 

T 

where is a stochastic variable that describes thermal noise. It has a 
symmetric Gaussian distribution with a two-point function 

(^M^W)} =2D{r)G e(l {T)8{T-T , )6{r,-r 1 l ). (2) 

The amplitude of the noise is given by G cci (t), which depends on the tem- 
perature and can be written as G cq (T) = q 2 Tn cq (T), where q is the elemen- 
tary charge and n eq (r) is the equilibrium particle density. 

The stochastic term in Eq. (1) can be eliminated by considering the 
two-point function G(t, r/ — r/') = (p(r, ??)/5(t, rf)). It satisfies the equation 
of motion 

d r G{r,) = ^^d 2 [G( V ) G cq (r)S( V )} . (3) 
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The Fourier modes G(k r) ) = J dr]e ik " v G(r]) satisfy 

2D(r)kl 

d T G(k n ) = [G(k v ) - G eq (r)] . (4) 

If we assume that initially at r = r; n i, the system is in equilibrium, i.e., 
G(Ti n i) = Geq(Ti n i), then the solution is 

G(r, k n ) = G cq (r) - f dr' e^^^C^r'), (5) 



where A 2 (r') = 4 /J", dt\D{r) /f 2 ]. 
Back in coordinate space, we have 

r p -1 2 /2A 2 (t') . 

G(r , 77) = G oq (r - dr> '^_==G cq (r') (6) 

Assuming that the charged particles are ultrarelativistic, their equilibrium 
density depends on the temperature as n eq oc T 3 . Entropy density scales in 
the same way, and therefore the temperature decreases with the expansion 
as T oc r -1 / 3 , meaning that G oq remains constant. This means that to 
a good approximation, G eq is simply a sum of delta functions at phase 
transitions and other non-adiabatic events. Consequently, G(r, 77) becomes 
a sum of Gaussians with different heights and widths. 

As an example, let us assume that the system is initially in thermal 
equilibrium in the QGP phase so that G eq (ri n i) = Gqgp- When the system 
enters the hadronic phase at r = T tr , G oq jumps discontinuously to Ghad, 
which is higher because the elementary charge q is ±1 instead of 2/3 or 
— 1/3. Thus, G cq (r) = (Ghad — Gqgp)S(t — r tr ), and the final two-point 
function is 

-r, 2 /2A 2 (r tr ) 

G(r, „) = G had %) - (G had - G QGP ) (7) 

^27rA^(T tr ) 

The charge fluctuation signal (AQ 2 ) can be written in terms of G(r, 77) as 

1 /-A77/2 

(AQ 2 ) = — / drjdr]' G{t, rj — rj'). (8) 

AJ 7 J-A v /2 

It is easy to see that it, indeed, gives the QGP result if Ar/ ^> A(r tr ) and 
the hadronic result otherwise. 

If we then assume that the hadrons become non-relativistic and start 
to annihilate at a later time r nr , the delta function peak spreads into a 
Gaussian 

-r, 2 /2A 2 (r„ r ) -r, 2 /2A 2 (r tr ) 
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It is interesting to note that at short enough distances (t]<A(T nl )), the 
correlator is positive. 6,7 . This also means that for small Ary, the charge 
fluctuation should behave as (AQ 2 ) oc Arj 2 . 

3. Long-range forces 

The discussion in Section 2 applies to global charges such as baryon number 
or strangeness, but for electric charges, one has to take into account the 
long-range Coulomb interaction, and the motion of the charges is not purely 
diffusive. By combining Ohm's law j = aE with Gauss's law V • E = p, the 
diffusion equation gets modified and becomes 

d T p = ^d 2 v p - <7(r)p + d^r,. (10) 
T 

It is instructive to note that when all the coefficients are time-independent, 
one finds that the equilibrium two-point function is 

which shows that the Debye screening mass mjj is given by m 2 D = a/D. 

Analogously to Eq. (3), we can derive an equation of motion for the 
two-point function 

d T G(k n ) = LlJL [ G (k n ) - G cq (r)] - 2a(T)G(k n ). (12) 

T 

Assuming that initially G(k v ) vanishes, the solution is 

G(r, k v ) = G oq (r) - f rfr ' e -^A 2 (r')^-2S(r) [ 2a ( r ') Geq ( r ') + Gcq ( T ') 
Jo L 

(13) 

where S(r') = J , dfa(f). Again, this is simply a superposition of Gaus- 
sians, and can therefore be easily Fourier transformed back to coordinate 
space. 

As an example, let us consider a simple case in which G eq jumps from 
to 1 at T tr = 1. Furthermore, we assume that D(t) = 0t, where j3 is 
constant, and cr(r) = tti 2 d D{t) with constant m 2 D . The correlator G(T 7 k v ) 
has been plotted in the left panel of Fig. 1 for m 2 D — [global charge, from 
Eq. (7)] and for m 2 D = 0.01 [local charge with long-range forces]. One can 
see that for a local charge, the correlator is more strongly peaked around 
zero. 

We then imagine that G eq drops instantaneously to zero at r nr = 10 
as the charged particles become non-relativistic. The correlators for global 
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Figure 1. G(r, k v ) as a function of k v at r = 10 (left) and at r = 13 (right). The two 
curves correspond to m? D = and m? D = 0.01. 

and local charges at r = 13 are shown in the right panel of Fig. 1. Because 
of the long-range forces, annihilation is faster and 

4. Conclusions 

We have seen that the charge correlators measured at late times carry de- 
tailed quantitative information about the properties of the system at differ- 
ent stages of the collision. Particle-antiparticle pairs produced later on by 
neutral resonances etc. will add an extra contribution to these correlations, 
but if it is properly understood it can be subtracted, at least in principle. 

The signals will be much stronger in global charges, i.e., baryon number 
and strangeness than in the electric charge, but unfortunately they are also 
much more difficult to measure. It is possible, though, that some of the 
interesting features survive even in the electric charge correlators, but a 
more detailed calculation is needed to find out if that is actually the case. 
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